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Abstract 



We consider biorthogonal polynomials that arise in the study of a generalization of two-matrix Hermitian 
models with two polynomial potentials Vi{x), V2{y) of any degree, with arbitrary complex coefficients. 
Finite consecutive subsequences of biorthogonal polynomials ("windows"), of lengths equal to the de- 
grees of the potentials Vi and V2, satisfy systems of ODE's with polynomial coefficients as well as PDE's 
(deformation equations) with respect to the coefficients of the potentials and recursion relations connect- 
ing consecutive windows. A compatible sequence of fundamental systems of solutions is constructed for 
these equations. The (Stokes) sectorial asymptotics of these fundamental systems are derived through 
saddle-point integration and the Riemann-Hilbert problem characterizing the differential equations is 
deduced. 
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1 Introduction 



In H] the differential systems satisfied by sequences of biorthogonal polynomials associated to 2- 
matrix models were studied, together with the deformations induced by changes in the coefficients of the 
potentials determining the orthogonality measure. For ensembles consisting of pairs of N x N hermitian 
matrices Mi and M2, the U{N) invariant probability measure is taken to be of the form: 

— dM(Mi, M2) — cxp Ur {-ViiMi) - ^2(^/2) + MiM2)dMidM2 . (1.0.1) 
Tn Tn n 

where dMidM2 is the standard Lebesgue measure for pairs of Hermitian matrices and the potentials Vi 
and V2 are chosen to be polynomials of degrees di + 1, d2 + 1 respectively, with real coefficients. The 
overall positive small parameter H in the exponential is taken of order when considering the large 
TV limit, but in the present context it will just play the role of Planck's constant in the string equation. 

Using the Harish-Chandra-Itzykson-Zuber formula, one can reduce the computation of the corre- 
sponding partition function to an integral over only the eigenvalues of the two matrices 



TN-= J J d^i{Mi,M2) cx / ]^da;,dy,A(a;)A(y)e-*^f=i^i(^^)+^^(^^)-^^''^ , (1.0.2) 

i—l 

and then express all spectral statistics in terms of the associated biorthogonal polynomials, in the same 
spirit as orthogonal polynomials are used in the spectral statistics of one-matrix models ||2^. In this 
context, what is meant by biorthogonal polynomials is a pair of sequences of monic polynomials 

7r„(a;) = x" + • • • , a^iy) = + ■ ■ ■ , nSN (1.0.3) 

which are mutually dual with respect to the associated coupled measure on the product space 

fdxdy ^„(x)a™(y)e-*(^^(")+^^(^)-"'') = /i„,S„„. (1.0.4) 



In this work, we use essentially the same definition of orthogonality, but extend it to the case of 
polynomials Vi and V2 with arbitrary (possibly complex) coefhcients, with the contours of integration 
no longer restricted to the real axis, but taken as curves in the complex plane starting and ending at 
00, chosen so that the integrals are convergent. The orthogonality relations determine the two families 
uniquely, if they exist ||l^, ^ . 

It was shown in ^, || that the finite consecutive subsequences of lengths ^2 + 1 and di -I- 1 respectively, 
within the sequences of dual quasi-polynomials: 

M^) = ^7r„(x)e-*^^(^) , My) = -^a„(y)e-*^^('') , (1.0.5) 

beginning (or ending) at the points n = N, satisfy compatible overdetermined systems of first order 
differential equations with polynomial coefhcients of degrees di and c?2, respectively, as well as recursion 



1 



relations for consecutive values of N. In fact, certain quadruples of Differential-Deformation-DifFerence 
equations (DDD for short) were derived for these "windows" , as well as for their Fourier Laplace trans- 
forms, in which the deformation parameters were taken to be the coefficients of the potentials Vi and V2. 
It was shown in |3| that these systems are Frobenius compatible and hence admit joint fundamental 
systems of solutions. 

In the present work we explicitly construct such fundamental systems in terms of certain integral 
transforms applied to the biorthogonal polynomials. The main purpose is to derive the Riemann-Hilbert 
problem characterizing the sectorial asymptotic behavior at a; = 00 or y = cx). 

The ultimate purpose of this analysis is to deduce in a rigorous way the double-scaling limits N 
00, hN — 0(1) of the partition function and spectral statistics, for which the corresponding large N 
asymptotics of the biorthogonal polynomials are required. (See |p^ , p2[ and references therein for further 
background on 2-matrix models, and ^, |lj, [l^, ^ ^ for other more recent developments.) The study 
of the large N limit of matrix integrals is of considerable interest in physics, since many physical systems 
having a large number of strongly correlated degrees of freedom (quantum chaos, mesoscopic conductors, 
. . .) share the statistical properties of the spectra of random matrices. Also, the large N expansion of a 
random matrix integral (if it exists) is expected to be the generating functional of discretized surfaces, and 
therefore random matrices provide a powerful tool for studying statistical physics on a random surface. 
(The 2-matrix model was first introduced in this context, as the Ising model on a random surface [p2[.) 

It has been understood for some time that the 1-matrix model is not general enough, since it cannot 
represent all models of statistical physics (e.g., it contains only the {p, 2) conformal minimal models). In 
order to recover the missing conformal models {{p, q) with p and q integers), it is necessary to introduce at 
least a two-matrix model The 1-matrix models are actually included in the 2-matrix models, since if 
one takes d2 = I, and integrates over the Gaussian matrix M2, one sees that the 1-matrix model follows, 
and hence may be seen as a particular case. 

It should also be mentioned that most of the results about the 2-matrix model (in particular those 
derived in the present work) can easily be extended to multi-matrix models without major modifications 
(see the appendix of (P). Indeed, the multi-matrix model is not expected to be very different from the 
2-matrix case (In particular, it contains the same conformal models.) 

The present paper is organized as follows: in Section]^, we present the required formalism for biorthog- 
onal polynomials, beginning with the systems of differential and recursion relations they satisfy, and 
recalling the main definitions and results of We then derive the fundamental systems of solutions 
to the overdetermined systems for the "windows" of biorthogonal polynomials in two ways: one, by ex- 
ploiting the recursion matrices Q,P for the biorthogonal polynomials, which satisfy the string equation, 
and another by giving explicit integral formulas for solutions and showing their independence when taken 
over a suitably defined homology basis of incquivalcnt integration paths. 
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In Section ^ we use saddle-point integration methods to deduce the asymptotic form of these fun- 
damental systems of solutions within the various Stokes sectors, and from these, to deduce the Stokes 
matrices and jump discontinuities at oo. The full formulation of the matrix Riemann-Hilbert problem 



characterizing these solutions is given in Theorem 3.1 



Acknowledgements: The authors would like to thank A. Kapaev for helpful comments concerning non- 



singular integral representations of the fundamental systems of solutions to eq. ( 2.1.3^ ) , which resulted 



in the addition of appendix A and the discussion at the end of Sec. 2.3 



2 Differential systems for biorthogonal polynomials and funda- 
mental solutions 

2.1 Biorthogonal polynomials, recursion relations and differential systems 

The definitions and notation to be used generally follow with some minor modifications, and will be 
recalled here. Let us fix two polynomials which will be referred to as the "potentials" , 

= uo + E ^^"^ ' - «o + E ■ ^^-^-^^ 

K=l .1=1 

In terms of these potentials we define a bimoment functional, i.e., a bilinear pairing between polynomials 
7r(x) and cr(y) by means of the following formula 



7r,cr 



dxdye-i'^^'^''^+^^^y^-''yK{x)a{y) . (2.1.2) 



The contours of integration r*^^\ F*^^) in the complex x and y planes remain to be specified. They will 
generally be chosen to begin and end at cxi, approaching it asymptotically in any direction that assures 
convergence. More generally, linear combinations of such integrals along various inequivalent contours 
may also be chosen. In fact there are precisely di (homologically) independent choices for the individual 
contours in the x plane and ^2 in the y-plane, due to the choice of the integrands |^ . The necessary and 
sufficient condition for the convergence of these integrals is that the contours approach oo in such a way 
that 

?StiVi{x)) +00, 3fi(V2(y)) +00. (2.1.3) 
For A; = 0, . . . 2di + 1 (mod 2di + 2) and / = 0, . . . 2d2 + 1 (mod 2^2 + 2), let us define the sectors 



r ^ ' ^'^^"^^ ^ 1,^^ + 2(rfi + 1) ' + 2(di + 1) 
■■= - arg(udi+i)/(di + 1) , 

■= - arg(wd2+i)/ (^2 + 1) • 
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Definition 2.1 (See Fig. 1) The wedge contours in the x-plane (resp. y-plane) (resp. ^j) for 
i = 0, . . . di (resp. j = 0, . . . 62), and the anti-wedge contours f (resp. f j^'^ for i = 0,...di (resp. 
j = 0, . . .d2) are defined as follows: 

• F-^^ (resp. ^'f^) comes from 00 within the sector S2i}_2 (resp. S2j^_2) and returns to infinity in the 

sector 62^^ (resp. S''2^). 

• fl"^^ (resp. f j^^J comes from 00 within the sector s'^^^ (resp. s'^^_i) and returns to infinity in the 

sector S'^l^ (resp. 

They have the property that 3ffVi(x) oo on the wedge contours T^^^ and ^Vi{x) —00 on the anti- 
wedge contours v'f'^ as \x\ 00 (and similarly for V2{y)). 




Figure 1: An example of wedge (solid) and 
anti-wedge (line-dot-line) contours in the x- 
planc for a case with di = 7 and real positive 
leading coefficient in Vi .The white sectors are 
the even numbered ones and the gray sectors 
the odd- numbered ones. (Only the first few 
are labelled). 



Note that since there are no singularities in the finite region of the x- or j/-planes for the integrands we 
are considering, we can deform the contours F^^-* and F^.^^* arbitrarily in the finite part of the x and y 
planes. Wc could therefore also take such contours as straight rays coming from the (2k — 2)-th sector 
combined with another ray going to infinity within the 2A;-th sector. Moreover, by the Cauchy theorem, 
we have the homological relation 

^ =0, ^'f^ = ' homologically . (2.1.6) 

i=0 j=0 
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There are therefore are only ^2 homologically independent contours F^^' and di of type F^.^^ (Similar 
remarks apply to the anti- wedge contours which will be needed later in this article.) 

In general, for any complex di x d2 matrix with elements {x*'-'} <=i,. .di define the homology class of 

j=l,...d2 

contours xF to be 



xF:=^x''^Ff' xFJ^^ (2.1.7) 



For brevity, we will denote the corresponding integral operator as follows. 



VV^^.^- / := / . (2.1.8) 



For a given x, we denote the corresponding bilinear pairing as: 

{Tr,a)^:= [ daidj/ 7r(a;)e-H(^i(^)+T^2(3')-^2')CT(y) , (2.1.9) 

and define two sequences of monic polynomials 7r„(a;), crn{y) of degree n such that they are biorthogonal 
with respect to this pairing 

(7r„,tTm)^ = hnSnm , (2.1.10) 
7Tn{x) = + ... , an{y) = 1/" + ... , (2.1.11) 

The matrix x*'^ must be chosen such that the N x N finite submatrices of the matrix of bimoments 

{Bij := {x\y^)}ij=o,...,N-i (2.1.12) 

are nonsingular for all N; i.e., the nondegeneracy condition that ensures the existence of the biorthogonal 
polynomials is given by 

A^(x):=det[(x',y^)J.^.^^ _^^0 ViV e N . (2.1.13) 

Since these A;v(>?)'s form a denumerable sequence of homogeneous polynomials in the coefficients x*'-' , 

this will hold on the complement of a denumerable set of hypersurfaces, and in this sense the permissible 

choices of x's arc "generic" , whatever the choice of potentials Vi and V2 ■ 

We introduce the corresponding normalized quasipolynomials and combine them into semi-infinite 

column vectors ( "wave- vectors" ) \l/(a;), $(y) as follows. 

00 00 

Mx) ■■= ^7r„(x)e-*^^(-) , My) ■■= ' 

^r(a;) := [iPoix), ...,Mx), -Y , $(y) := [My), ■.■,My), -Y ■ (2-1-14) 
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In matrix notation the biorthogonality reads 



dxdye " \p(a;)<i) (y) = 1 , 

.p oo oo 

where 1 denotes the semi-infinite unit matrix. 

We denote the Fourier-Laplace transforms of the wave- vectors along the wedge contours by 

¥-'\y):^ [ dxe^^^/\x) , i = 0,...,di 



(2.1.15) 



dye'' $*(?/) , j = 0, ...,d2, 



(2.1.16) 
(2.1.17) 



viewed as row vectors, with components denoted ip^^Hu) a-nd 0|^"'-'(a;) respectively. Because of eq. (2.1.6) 
only di (or ^2) are linearly independent 



j=0 °° 



(2.1.18) 



1=0 



The recursion relations for the biorthogonal polynomials are expressed by the following matrix equa- 
tions for the wave- vectors [0: 



X vl/(a;) = Q ■^{x) , hdx '^{x) = —P vl/(x) , 

oc 00 00 00 

y$(y) = P*$(y) , hdy <^>{y) = -Q' ^{y) , 



(2.1.19) 
(2.1.20) 



where the matrices Q and P have finite band sizes ^2 + 1 and c?i -f 1 respectively, with Q having only 
one nonzero diagonal above the principal one and P only one below 



Q 



ao(0) 7(0) 
ai(l) ao(l) 7(1) 








P 



ad2(c^2) ad2-i(c^2) 


/?o(0) /3i(l) ••• PM) 

7(0) Ml) /3i(2) ■■. PdAdi + 1) 

7(1) /3o(2) 

7(2) /3o(3) ■■• 



(2.1.21) 



(2.1.22) 



These semi-infinite matrices satisfy the string equation 

[P,Q] = hi 



(2.1.23) 
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(See 1^, |T^, Q for simple proofs of these assertions.) For the dual sequences of Fourier-Laplace transforms, 
a simple integration by parts in eqs. ( ^.1.19| ), ( 2.1.20[ ) gives 



X (x) = {x)Q , hd^ {x) = {x)P , j = {),..., d2 

OO OO OO OO 



(2.1.24) 
(2.1.25) 



Notice that integration by parts is allowed due to the exponential decay of the integrand along the chosen 
contours. 

We recall from |^ the definition of dual sequences of windows. 

Definition 2.2 Call a window of size di + 1 or d2 + 1 any subset of di + 1 or d2 + I consecutive 
elements of type ipn, (/) j 4'n or ip , with the notations 



^! ■.= [lpN^2^---,'lpN]\ N>d2, $ := [(/)Ar_rf,,...,0Ar]* , N > di 

N N 

N f N 

^ := [-0Ar_l, . . . , V'W-kii-l] , > 0, $ := [0Ar_i, . . . , (l)N^2-l] : N>1 



N 



'In 



N>d2, , N>d, 



N 



(2.1.26) 

(2.1.27) 
(2.1.28) 



> 1, 



N+d2-V ' 



N > 1 .(2.1.29) 



Note the difference in positioning and size of the windows in the various cases, and the fact that the 
barred quantities are defined as row vectors while the unbarred ones are column vectors. The notation 

N N+1 

we are using here differs slightly from that used in H], in which, e.g., $ here would be denoted $ in B. 

N N 

Definition 2.3 Each of the two pairs of windows (\p, $) and (<I), ^) of dimensions d2 + 1 and di + 1 

N N 

respectively, will be called dual windows (and these are defined for N > d2 and N > di respectively) . 

We now recall the notion of folding, as defined in Q . We give here only the main statements without 
proofs and refer the reader to ref. for further details. Let us introduce the sequence of companion-like 

N 

matrices a(x) and a(x) of size do + 1 

N ~ 



a.(x) 

N 



a{x) 








j(N) 
x~ac(N) 
7(JV-1) 

-ai(jV+l) 
7(JV-1) 



-ad^(jV+rf2) 
7(iV-l) 





1 

-ai{N) {x-ao{N)) 
fiN) 7(W) 

10 0" 

•• 

1 





TV > do 



(2.1.30) 



N>1, 



(2.1.31) 
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N 

and also the analogous sequence of matrices h{y) and b(y) of size di + 1 

N 



1. 

N 







1 














-/3di(iV) 



N 

hiy) 





1 

-Pi(N) iy-f3o{N)) 

1 



N>di 



•• 

1 



, N>1. 





The first equations in (|2. 1.191) , ( |2.1.20|) and (|2.1.24|) , ( ^.1.25|) imply the following 

N N 

Lemma 2.1 The sequences of matrices a, b and a, b implement the shift N N + 1 and N i 

N N — 

in the windows of quasi-polynomials and Fourier-Laplace transforms in the sense that 



a ^{x) = V]/ (x) 

N N N+1 

h^iv) = $ (y) 

N N N+1 



N N+1 jv 

$(a;) = ^ (x) a 

N N+1 N 

My) = * (y) b 



and in general 



(2.1.32) 



(2.1.33) 



N-1 

(2.1.34) 
(2.1.35) 









N 




~ a 


• a^f , 


$ 


N+i 


N+j-l 


N N 


N 


$ 


= b • 


•b$ , 




N+i 


N+j-l 


N N 





N 



N 

h 



(2.1.36) 
(2.1.37) 



N N 



where ^ , $ here denotes a window in any of the Fourier- Laplace transforms defined in eqs. (2.1. It ), 



(2.1.17). 



(In this lemma, and in what follows, if no superscript distinguishing the integration path in the Fourier- 

N N 

Laplace transform in eqs. ( 2.1.16| ), (2.1.17) is present in $, this means that the result being discussed 

holds for all cases.) 

Equations (|2.1.34[) , ( ^.1.35| ) will henceforth be referred to as the "ladder relations" . We refer to the 
process of expressing any ipn (x) or (j)^ (x) by means of linear combinations of elements in a specific window 
with polynomial coefficients, as folding onto the specified window. We also recall that matrices ^(y) and 

N N J-, 

a{y), h{y) and h{y) are invertible (see |g]). The opposite shifts are therefore implemented by the inverse 
matrices and the folding may take place in either direction with respect to polynomial degrees. 
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N N 



Lemma 2.2 (From with the adapted notation.) The windows of quasi-polynomials \I>(a;), $(?/) and 

N N 

Fourier-Laplace transforms ^{x),^{y) satisfy the following differential systems 

N>d2 + l, 



where 



^9 , , 
ax N 



N 

Di{x)'qj{x) 



d N 



^9 , , 
oy N 
d ^ 



N 

N N 

^ix)D^(x) 



N 

D2{y)^{v) 

N 

N N 



di 



N N N ^ N 

Di{x) :=/3_i a -i+/?o+V/3, a 

N-l ^ ■' N+j-l 



J = l 



N>d2 + l 

N>di + l 
N>di + l 



a • • • a e gld2+i N 

N+j-2 N 



N 



N 



N 



N 



di 



i=i 



N-lN-2 N-] N , r 1 

a a • • • a /3. e gld^+iM , 



AT 



diag [P,{N+]-d2),P,{N+j-d2 + l), . . • ,/3,(iV+j)] , j = 0, . . . di 



N 



P_^:= diag [-i{N-d2-l), ...,l{N~ 1)] 



N 



B:=Ai&g{Pj{N -l),Pj{N),...p,{N+d2~l)) , j = 0, , 



N 



:= diag (7(iV - 1), 7(A^), • • ■ , 7(^ + ^2-!)) 



AT 

^ , S ^ -1 ^ \ ^ 

D2(y) :=a_i _b +ao+2_^ 



N-l 



N , r 1 

aj b b •••be gldi+i[y\ ■ 

N+j-1 N+j-2 N 



N 



d2 



M.2[y)-^h a_i+ao + 2^ b b ••• b a^- e gldi+i[y\ , 



3 = 1 



N 



N 



diag[aj{N+j-di),aj{N+j-di + l),...,aj{N+j)] , j = 0, , 



a_i:= diag[7(iV-di-l), . . . n{N - 1)] 

N 

aj:=dieig{aj{N~l),aj{N),...aj{N + di-l)) , j = 0, . . . , da • 
«_i :=diag(7(iV-l),7(iV),...,7(7V+di-l)) . 

Moreover the systems ( 2.1.3^ )-I ^.1.41 ) are compatible with the ladder relations (2.1.34), (2.1.31 



N+l N 

Di (x) = a.(x) Di(x)sl^'^(x) — h a'(a;)a~-^(a:;) 

AT N N N 

N+l 



N 



N 



N 



N,, , N_ 



Di (x) — s.{x) Di{x) a {x)~ha'{x)a (x) 
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N+1 JV 

D2 (y) = h{y) D^{y)h-\y) -hh'{y)h-\y) 

N JV N N 

N+1 N N N , JV JV , 

D, {y) = h{y)D,iy)h-\y) - hh'{y)h-\y) 



(2.1.56) 
(2.1.57) 



Remark 2.1 We stress that the windows taken from any wave-vector solution to eqs. 

and eqs. (2.1.24), ( 2.1.2^ will automatically satisfy the systems specified in Lemma 2.'i (differential and 

difference equations). 

In the following, wc will only consider the fundamental systems of solutions to the sequence of equa- 
tions (2.1.3S), (2.1.39), (2.1.54), (2.1.55), since the corresponding dual sequence of equations for the 

JV 

c^uantities $(?/), ^{y) may be treated analogously by just making the appropriate interchange of no- 
tations X 
planes. 
In 



JV 



y, ^ di d2, etc., and the interchange of integration contours in the x and y 



it was shown that there exists a natural nondegenerate Christofjel-Darboux pairing between 
~ _ JV 

any pair of solutions '^{x) and $(a;) to eqs. (2.1.34) and (2.1.39) such that 



^jv _^ 



JV 



JV . 



JV 



$(a;) A*(a;) 

~ JV 



JV 



(2.1.58) 



is constant both in x and N, where the invertible matrix A defining the pairing is 

















^(7V_1)1 


JV 






a2iN) 







A := 









a2{N + l) 































adiN+d2-l) 






(2.1.59) 



JV N 

This follows from the fact that the matrices Di and are conjugate to each other by means of the 

JV 



matrix A (Theorem 4.1 in 



JV AJ N N 

KDi{x) = D^ix)A , 



(2.1.60) 



JV 



together with the shift relation ( 2.1.54 ). The {d2 + 1) x (^2 + 1) matrix A is the only nonzero block in the 



commutator 



JV-l 

n ,Q 




JV-l 



, where n denotes the projector onto the first N basis elements (a "canonical" 



JV 



projector). In the following it is convenient to use the same notation A both for the finite matrix and 
the semi-infinite one. 

_ JV 

It was also proved in ||] (Proposition 3.3) that one can choose the windows ^ and $ as joint solutions 
of the PDE's following from the infinitesimal changes of the coefficients {uj} and {vk} of the potentials 
(deformation equations). For this choice, the pairing also becomes independent of these deformation 
parameters. 
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2.2 Fundamental solutions of the Di and D^^ systems 



N 

In this section we explicitly construct solutions of the pairs of dual ODE's defined by the matrices Di{x) 

N 

and D.i{x). In fact these solutions will simultaneously satisfy the deformation equations and the ladder 
recursion relations in A^. It will be left to the reader to formulate the corresponding statements for the 

N N 

other pairs {D2{y) ^ Djiiv)) i which are essentially the same, mutatis mutandis. 



As mentioned in Remark 2.1, taking windows within any wave- vector solution to eqs. ( 2.1.19| ), (2.1.24) 



one obtains solutions to the difference-differential equations in Lemma 2.2. Therefore we could try to 
construct ^2 + 1 such wave- vector solutions in order to obtain a fundamental system for eqs. ( 2.1.39| ). 



On the other hand it will be shown in Proposition 2.1 below that the relations (2.1.24) have precisely 3,2 
linearly independent wave- vector solutions (given by the d2 Fourier-Laplace transforms ( ^.1.17| )) while 



relations (2.1.19) have only one solution (given by the quasipolynomials) . Therefore we will look for 
solutions to slightly modified recursion relations differing from ( ^.1.19| ), ( |2.1.25 ) by changing the initial 
term^. The solutions to these modified recursion relations {'0n(a;)}n>o or {V'„(2;)}n>o will still satisfy 
the unmodified recursions relations, but only for n large enough. This will provide us with the required 



solutions of the difference-differential equations of Lemma 2.2 



We start by proving that the Fourier-Laplace transforms and the quasipolynomials are the only 
wave- vector solutions to eqs. ( 2.1.19| ), (2.1.24). 

Notational remark. Throughout the remainder of this section, ^!{x) and ^{x) will be used generally 
to denote arbitrary wave-vectors consisting of arbitrary solutions and not just the previously defined 
quasipolynomials and their Fourier-Laplace transforms. 



Proposition 2.1 The semi-infinite systems 



X <i!{x) = Q \l/(x) 

oo oo 

hdx ^{x) — —P vl/(x) 



(2.2.1) 



X $(x) = ^{x)Q 

oo oo (0 0\ 

hdx^{x)^^{x)P ^ ' 

oo oo 

have 1 and o?2 linearly independent solutions, respectively, for the wave-vectors '^oo{x) o,nd ^^{x). (A 



similar statement holds for $, as solutions to eqs. (2.1.2L ) and (2.1.25) with d2 di). 

oo oo 



Proof. The compatibility is guaranteed by the string equation (2.1.23). Recalling [|| that 



p - vm 



>0 



Q - vi{P) ] = 

<o 



(2.2.3) 



^This is not new in the context of^ttbogonal polynomials, where there exist, besides the orthogonal polynomials, 
solutions of the second kind. (See e.g. [pi Ea]). 
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we have 



= 







(P - VliQ)) VI/ 

OO 



(2.2.4) 
(2.2.5) 



Thus we have only one solution of the first equation (up to normalization) and d2 independent solutions 
for the second. Using the x recursion relations for the V' sequence, we can build the rest of the sequence by 
starting from the first term -00 (2;) , since the matrix Q has only one nonzero diagonal above the principal 
one. On the other hand, using the dx recursion relations we can build the rest of the sequence starting 



from any given solutions (p^ of (2.2.5) because the matrix P has only one diagonal below the main one. 
Q. E. D. 



The solutions of Propostion 2.1 are, up to multiplicative constants, exactly the solutions given by 
the quasipolynomials for the wave-vector v[/ and the different possible $'s corresponding to the different 



solutions of eq. (2.2.5), which are just the ^2 different Fourier-Laplace transforms of the quasipolynomials 



i(y). Indeed, the solutions 4>Jyx) of eq. ( 2.2.5 ) can be expressed by 



dye 



k(^v-v.,(v)) 



(2.2.6) 



iv) 

where F is any of the contours or a linear combination of them. 

As indicated in the introduction to this section, in order to find a fundamental system of solutions to 



the Difference-Differential equations in Lemma 2/1 , we consider systems for the wave- vectors modified by 
the addition of terms that only change the equations involving the first few entries of the wave-vectors. 
The recursion relations will therefore remain unchanged for n sufficiently large. 



Proposition 2.2 The semi-infinite systems: 

{ x-^{x)=Q-^{x)-W2{-hdx)F{x) 

] h^,^^,{x) ^-p^^,{x) + F{x) 



x^{x) ^ ^{x)Q + U{x) 
hdx^ix)'^ ^{x)P + U{x)Wi{x) ' 

00 00 



(2.2.7) 



(2.2.E 



where 



P - hdx 



(2.2.9) 
(2.2.10) 



W.ix):=Yi^-^ 

Q - X ^ 

F{x) [/(a;), 0, 0, • • •]* , U{x) := [w(x), 0, 0, • • •] , 

toth ha.e + 1 .nearly independent solutions for the extended unknown wa.e-.eetors f /(.), .,(.)) and 

\ 00 

( u{x),^{x 
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Remark 2.2 The terms in the RHS of eqs. i\2.2.1[ ), ^.2.^ which are not present in eqs. ^2.2.:^ , i\2.2.^ 
are semi-infinite vectors with at most d2 nonzero entries for ( 2.2.'\ ) or di for (2.2.1 ). Indeed the matrices 
Wi and W2 (which are polynomials in Q and P of degrees rfi — 1 and ^2 — 1 respectively) have finite hand 
sizes, and hence when acting on the vectors U{x) and F{x), which have only a first nonzero entry, produce 
a finite number of nonzero entries. As a result, windows taken from any solution of the modified wave- 
vector systems ( 2.2.7 ), ( 2.2.1 ) will provide solutions of the Difference- Differential systems of Lemma 2.i , 
but only beyond a minimal N value (di or ^2 respectively) . 



Proof. The compatibility of these systems is not obvious. We have: 



(2.2.11) 
(2.2.12) 
(2.2.13) 
(2.2.14) 



Since only the first entry of F is nonzero and Q — V^iP) is strictly upper-triangular the second term 
vanishes. The last term gives the following ODE for the first entry of F{x) 



^Q[-P^+F] -d,W2i~nd,)F + P [Q-if^W2{-hd^)F] +xF 
= [P, Q] * +QF - [hd. + P)^l(^l_i|Iz^F(x) + xF 
= /i* + (0-l/2(^)V+ (x-V^i-hd^UF. 



vi{-hd,)f{x) - xf{x) 



(2.2.15) 



The solutions of cq. (2.2.15) are easily written as Fourier-Laplace integrals, and give a compatible system 
(see eq. (2.3.16), where we also fix the most convenient normalization for them). We denote them by Z*^"-* 
with a = 0, . . . , c?2, with /'■"^^ = the trivial solution, corresponding to the unmodified system (2.2.1). 
Fixing any such solution /(x), we can now solve for First of all, one can prove by induction that 



^^,^x^^ +^ — —W2i-hd.^)F. 



(2.2.16) 



Next we compute as for the previous proposition 



= 



(P^VliQ))^ 



Moi^) - v{ix)Mx) + [1 - Wi{x)W2{~hd,)U fix) 



(2.2.17) 
(2.2.18) 



Thus, ^^0(2;) must solve this first order system of ODE's. Since there are d2 + I choices for the function 
/ = f'""\ we correspondingly obtain d2 + 1 independent solutions 'I'("^(a;) to the system. 
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Now consider the second system ( 2.2.^ ). The compatibihty gives 

h^ = [hd^,x]^^hd^ (^Q + u) - x(^P + UWi{x)^ (2.2.19) 

^P + UWi{x) ] Q + hd^U- ( $Q + C/ ) P-UWi{x)x (2.2.20) 

(2.2.21) 

(2.2.22) 

Since the first entry of U is nonzero and P ~ V({Q) is strictly lower-triangular, the second term vanishes. 
The last term gives the following ODE for the first entry of U (x) : 

hu'{x) = Vl{x)u{x) =^ u{x) = cei ^^^""^ or u{x) = . (2.2.23) 

We next consider the solutions By a computation similar to the previous one, we have 



m Q]-UP + c/-llM-iiM(g -x) + hd.u 

oo ^ X 

h^^ufp - vliQ)] + (hd^ - v{{x)]u{x) 



x^-U -V^ihda,)^ +UWi {x) 



X - V^ihd^) ]Mx) - uix) 1 - Wiix)W2ifid^) 



vi{P)-vi{hd, 

P-fidx 



(2.2.24) 



(2.2.25) 



This is a (i2-th order inhomogeneous ODE for the function (j)^[x). Choosing u{x) = as solution to 
(2.2.25) we have 3,2 independent solutions corresponding to the "unmodified" system ( |2.2.2 ). The choice 

u{x) = ceT^Vi{x) (2.2.26) 

leads to one more independent solution of the inhomogeneous equation. Q. E. D. 

We denote the solutions to system ( 2.2.8 ) in general as a = 0, . . . ^2, where a = corresponds 

oo 

to the inhomogeneous solutions and a = 1, . . . (i2 correspond to the independent solutions to the homo- 
geneous case {U{x) = 0). We will give explicit integral representations for the (i2 + 1 solutions in the 
next section. With the solutions ^p'-"-', ^^^\ a, /3 = 0, . . . , c?2, we can construct the (^2 + 1) x {d2 + 1) 

CO oo 

modified kernels 



Af-l 



k[^/\x,x') y 4>^''\x)^\^\x') = ^^''\x)W(P\x') 

— " — oo 



n=0 



(2.2.27) 



and also obtain the following modified Christoffel-Darboux formulae (cf. 
Proposition 2.3 (Christoffel-Darboux Kernels) 

{x-x')Kt/\x,x') 



(2.2.28) 



Af-l 



t(")(a^)Q + '5aoC/(^) n ^^P\x') - ^^''\x) n [Q^>^^\x')~W2(~hd',)F^'^\x')] (2.2.29) 



oo 



N-l 

n ,Q 





*('')(x') + -JaO"(a;)Vr(2;') +t^"Ha^)'nV2(-?i5i)^^('^)(a;') . (2.2.30) 
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N-1 



Recall that the commutator of the finite band matrix Q with the projector n gives a finite-rank semi- 



N 



infinite matrix that corresponds to the ChristofFel-Darboux kernel matrix A- As a corollary, we obtain 
the pairing between the solutions of the systems ( 2.1.39| ) by setting x = x' 



Corollary 2.1 



N 



N 



N 



N 



N 



N-1 

n ,( 



N-1 



= S„ou{x)4^\x) + ^^"\x) n W2hhd,)F^f'\x) 



(2.2.31) 



As shown in this is a constant (in x) if N > d2. But if N > d2 the projector in the second term is 
irrelevant because the vector W2{—hdx)F^^\x) has only its first d2 + 1 entries nonvanishing, and hence 
we have, for N > d2, 



N J 00 P + ndx 

If a 7^ then 

$(")(a;)p = hdx^'^"\x) , 

00 00 

and hence we have, for N > a 7^ 0, 



00 00 00 hd X + fid X 



(2.2.32) 



(2.2.33) 



hdx' + hdx 



In this expression i/i^ and / are kernel solutions of the pair of adjoint differential equations 



{v^ihdx) - x)ci,l^\x) = , {v^i-hdx) - x)f^^\x) = , 



(2.2.34) 



(2.2.35) 



and the last expression in eq. (2.2.34) is just the bilinear concomitant of the pair (which is a constant in 
our case). For P — we have 



$(")(^) Avl/W(a;) = 5^ou{x)^l°\x) = 5„oce*^i(-)-^e-*^i(-) = 6^0^ 

00 00 V'iO V"0 



(2.2.36) 



For completeness we recall that the matrices P, Q satisfy the following deformation equations 



Q,U 



K 



K 



hdyjQ = 
hdyjP = 



P,V"' 



(2.2.37) 
(2.2.38) 
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where 



U 



-^{iQ^'ho + lmo} ' ^'■■=~-j{[P']<o + ln} • (2-2.39) 

(Note that in the matrix that here is denoted by P was denoted —P.) Here the subscripts < (resp. 
> 0) mean the part of the matrix below (resp., above) the principal diagonal and the subscript denotes 
the diagonal. 



Proposition 2.4 The equations in Proposition 2.i are compatible with the deformation equations 



hdu^F{x)^\j'^F{x) 



oo 

hd,,F{x) = 

OO 

hduit U {X) : 



-v-^^Gx) - 
J 

- $(x)U^ 



1 p' 



J P + hd^ 



Fix) 



1 .Q^'-x^ 



K 



U 



Q^x 

U{x) 



(2.2.40) 



(2.2.41) 



(2.2.42) 



(2.2.43) 



[ hd^, $(x) = ^{x)\-^ 

< OO OO 

[ hdy,U{x) = U{x)V^ . 

Proof. We only sketch the proof, which is quite straightforward but rather long. 

Compatibility of the deformation equations amongst themselves follows from the zero curvature equa- 
tions 



(2.2.44) 



(which are established in the general theory |^), together with the fact that both vectors U{x) and 
F{x) have, by assumption, nonzero entries only in the first position. 

Compatibility with eqs. ( ^.2.7 ), (2.2.8) requires some more computation. We first prove the com- 
patibility between eq. (2.2.8) and the two eqs. (2.2.42), (2.2.43). Compatibility between the equations 
involving multiplication by x and application of dyj'. 



hxdyj ^(x) = $(x)QV"' + U{x)\-^ 

OO OO 

hd^,x^{x)=^{x)Y'^Q + ^{x)[Q,V'] + U{x)Y^ = RHS of (2.2.45) . 

cx: OO OO 

Compatibility between the equation involving multiplication by x and application of : 

1 qK _ K 

xhdu^ ^{x) = ~ $(x)QU^ - U{x)\]^ + —U{x)^^—- X 

OO OO JV Ccj/ - X 



(2.2.45) 
(2.2.46) 

(2.2.47) 
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Q - X 



$(x)U^g + — 



Q-$(a;)[Q,U^]+ — +U^o C/(a;) 



K 



(2.2.48) 
(2.2.49) 



1 „^ 1 ,Q^-.T^ 



$(x)QU^ + -;7(x)(Q^-a;^) + -C/(x) 
oo A A Q - a; 



a; + ( ^ + U^o ) U{x) (2.2.50) 



(*) 



RHS of ( p^2l7| ), 



(2.2.51) 



where in the step (★), we have used the fact that (smce only the first entry of U{x) is nonzero) 



-U{x)Q'' + ^^oU{x) = -U{x)\]^ . 
Compatibihty between the dx and du^ equations is a bit more involved: 
h^dAj $(x) = hdx ^{x)Y'' = $(x)PV^ + U{x)Wi{x)W' 

OO OO OO 

h^d^jdx^ix) = hdyj ( $(x)P + f/(x)VKi(x) ) = 



= $(a;) V'P + $(a;) [P, V"'] + U{x)Y-^Wi (x) + U{x) [Wi {x) , V-^] 

oo oo 

Finally we have the dx^ dyj compatibility: 



(2.2.52) 



(2.2.53) 



RHS of (2.2.53) (2.2.54) 



h'dxd.,, t(x) = fidx { - $(a:)U^ + -U{x) 

OO 



Q - X 



^{x)P + U{x)Wi{x) + —hdxU{x) 

oo J K 

^{x)P + U{x)Wi{x)] + ^{hdx - Vi{x))U{x) 

oo J K 



l.U{x)hdx^-^ 



Q - X 



1 / lower tri. 

+ (WW - WW) + W(Q) - P 

+—U{x)P- 



Q^-x^ 
Q - X 

~ x^ 



X 



Q-x 



(**) 



^{x)P + U{x)Wi{x)] - l-U{x)Wi{x)iQ'' - x^) + l[/(a;)^^^P 
oo J K K Q - X 

+^{hdx - v({x))u{x)^-^^ + ^u{x)ihdx + hdg) ^ _ ^ 



K 



$(x)P - U[x)Wi{x) ( j^Q^ + 



1 



jUixW.ix)^^Uix)^ 
+ ^if^dx - y/(x)){/(x)^^^ + hUix) ^ ^ 



■P 



1 



$(x)P - C/(a;) ( + ) M^i(x) + '^U{x)Wi{x) + ^[/(x) 



1 



-P 



fiJ7(a;)- 



t/(a;)[Tyi(a;),U^] 
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K 



u, 



1 g^-:E^ 



+ 17(ft5,-y/(x))C/(a;) 

A 



C/(x)VKi(x) + -C/(a;) 

A 

;iC/(a;) 



-P 



U{x)[Wi{x),\5'^] . (2.2.55) 



_ — X 1^ — X 

In step we have used the fact that the commutator of P with a function of Q is equivalent to the 



(formal) derivative hdq^ while in step (-k-k*) we have used eq. (2.2.52). Computing the derivatives in the 
opposite order we obtain 



h^d^,d^^{x) = hdyj ^{x)P + U{x)Wi{x) 

oo \oo / 

1 nK_^K\ /^K 



-(-tWu« + ^c/(.)-^_ 



-hu{x[ 



X 



X 



U{x)[W,ix),\J 



X 

K 

K 



(2.2.56) 



Subtracting eq. (2.2.55) from eq. (2.2.56), we are left with 



= RHS of (2.2.55) - RHS of ( 2.2.5e ) = —{hda, - Vl{x))U(x 

K 



Q 



K 



{hd^~vi(x))u{x)^o , 



(2.2.57) 



where the implication follows from the fact that only the first entry of U{x) is nonzero. This is precisely 



the same compatibility condition that we found in Prop. 2.2. The other compatibility checks are also 
rather long but routine and are left to the reader to fill in. Q. E. D. 

2.3 Explicit Integral representations for the dual wave-vectors 

Proposition 2.5 The + 1 semi-infinite wave-vectors and functions u*^"-' (a;), a = 0, . . . ^2 defined 

oo 

hy 



p-j-{Vi(s)-sy) 

$(0)(a;) =e*^i(^) dsdy- $*(?/), u'-^^x) 



dye"^ $*(?/) , u^''\x) = 0, 



■,d2, 



(2.3.1) 



(2.3.2) 



are independent solutions of the modified system (2.2.t ). 

Proof. There is nothing to prove for ^^'^'j since for fc > these are just the Fourier-Laplace transforms, 



which satisfy the corresponding unmodified system (2.2.1), which is the same as the modified system 
with u{x) = 0. Let us therefore consider We first check the x recurrence relations. 



a;$(°)(a;) =e*^i(^) / dsdy 



(a;-s)e-*(^i(^)-^^) , e"*(^i(^)-^^) ' 

$ (y) + s$ (y) 
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jdsdy hdy<i>\y) 



OO 

where in (★) we have used eq. ( ^.1.201 ) and defined the semi-infinite vector 

t/(0)(a;) := [y^e^^^^^), 0, 



It is easy to prove by induction in k, starting from the first part of eq. (2.2.8), that 



(2.3.3) 



(2.3.4) 



We now consider the dx differential recursion. After shifting the derivative in x to one in s inside the 



integral, and integrating by parts, using eq. ( 2.3.4 ), we get 



(2.3.5) 



This proves that the given integral expression is indeed the additional solution to eq. (2.2.8). Q. E. D. 

Remark 2.3 We would also have solutions for any other choice of admissible contours (or linear com- 
binations of them) in eq. {2.3.i). This arbitrariness will be used later. 

Remark 2.4 The functions (ji'^^x) are piecewise analytic functions in each connected component of 
\ Uj=i Tj^'', but can be analytically continued from each such connected component to entire functions 
by deforming the contours of integration. 



With regard to the deformation equations studied in Propostion |2.4| , we have: 

Lemma 2.3 The wave-vectors ^/''^'^{x) (i.e. the quasipolynomials) and^^^\x), k — l,...,d2 (i.e. the 

oo oo 

Fourier-Laplace transform of the quasipolynomials 4>n{y)) satisfy the following deformation equations 

hdu^^'^''\x) = U^*(°)(a;) , (2.3.6) 

OO OO 

hd,,m^'^\x) = -V^*(°)(2^) , (2.3.7) 

OO oo 

hdu^^^^\x) = -$('=)(a;)U^ , (2.3.8) 

OO OO 

hd,,,^^^\x) = ^^^\x)Y-\ (2.3.9) 



Proof. This follows straightforwardly from the two equations (2.2.37), (2.2.38) together with the fact 
that taking a Fourier-Laplace transform commutes with differentiation with respect to the deformation 
parameters Q. E. D. 

We also remark that the given integral expressions satisfy the modified deformation equations given 
by the following proposition. 
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Proposition 2.6 The dual wave-vectors defined in Proposition 2.t satisfy the following deformation 
equations 



OO oo 



(2.3.10) 
(2.3.11) 



oo oo 

where the semi-infinite vector U^'^\x) is defined in eq. 

Proof. As for the {wj} deformation equations for tlie proof follows from Lemma pTs] and from the 

oo 

fact that the integral transform defining this wave- vector does not depend on the wj's. On the other 



hand, the uk deformations give (again using Lemma 2.2): 



hdu^ (x) = - $ (°) (a;)U^ + - $ (x) (x^ - Q^) 

OO oo K oo 



(2.3.12) 



from which the proof follows using eq. (2.3.4). Q. E. D 



Remark 2.5 Note that from the definition of the vector U{x), one can verify directly that it satisfies the 
deformation equations (2.2.4^), ( 2.2.4-i\ ) by using the equations for the normalization constant hg: 



ho — j dxdye 

h 

2 



-\{Vi{x) + V2{v)-xy) 



duK ln(/io) = , ^d.^j ln(/io) = 



(2.3.13) 
(2.3.14) 



Remark 2.6 In particular, the sequences of functions defined in Proposition 2J^ satisfy identical defor- 
mation equations for N large enough (N > di) and hence the relevant windows satisfy the full DDD 
equations specified in ^ for the barred (dual) quantities. 

We conclude with a discussion regarding the construction of an integral representation for the wave- 



vectors V]/. For the unmodified system (2.2.1) the solution is explicitly given by the quasipolynomials 

oo 

^I'^^^x). The other solutions of the wave-equations (2.2.7) are in 1 — 1 correspondence with the d2 



solutions of eq. ( ^.2.15| ) via eq. (2.2.18). We could try to define an integral representation of the form 

■^(z) , (2.3.15) 



corresponding to the following solution of eq. (2.2.15) 



2i'Kh Jfi 



t-y 



k — 1, . . . ,d2 



(2.3.16) 



However formula (2.3.15) is not well-defined, since the inner integral defines only a piecewise analytic 



function with jump discontinuities along the L^ s, and each of the anti-wedge contours crosses 
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at least one of the contours T^-^''. One could alternatively take the analytic continuation of the inner 
double integral from one connected component of Cj, \ IJ^Li T^,^'' to an entire function, which amounts to 
deforming the contours in xT to avoid any jump discontinuities. But then the outer integral is in general 
divergent, due to the behaviour of the integrand as y — > oo. 



Thus, as it stands, eq. (2.3.15) may only be viewed as a formal expression which at best defines a 
divergent integral. However, treating it naively (as though it were uniformly convergent with respect to 
the parameters {uj,vk}, as well as x) we would find that it indeed satisfies eqs. ( ^.2.71 ) (|2. 2.401) and 
( ^.2.4l[ ). Nonetheless we know that there are solutions vp^*^^ of eq. (2.2.7) ( 2.2.40| ) corresponding to the 



f'^'^^x) defined in (2.3.16), determined recursively from the solutions of eq. (2.2.18). Such solutions are 



uniquely defined by the choice of a particular solution to eq. ( 2.2.15 ), up to the addition of a solution of 
the homogeneous equation (that is, the quasipolynomials) forming ^!^'^\ 

oo 

There does exist a way to arrive at convergent integral representations, however, by forming suitable 
linear combinations of Fourier-Laplace transforms of a piecewise analytic function defined along appropri- 
ate contours. This was done for the case of cubic potentials by Kapaev in |2l[. In Appendix A, it will be 
indicated how the procedure used there may be extended to polynomial potentials of higher degree. For 



present purposes, however, we may view the solutions as just defined by the recursion relations (2.2.7), 



together with eq. ( |2.2.18| ). It follows from compatibility that these satisfy the deformation equations as 
well. 

As a further remark, we note that the solutions vj/^'^' ( modC vf^^-*) associated to the /''^^s defined by 



eq. (2.3.16), may be associated to any linear combination of the anti-wedge contours, used as contours 
of integration in ( ^.3.16 ) . This freedom will be used in the next section. 



2.4 Diagonalization of the Christoffel— Darboux pairing 



We know from Proposition 2.t that the windows constructed from the integral representations in Propo- 



sition 2.5 provide solutions to the DDD equations for the barred system (for > di). On the other 
hand, we are guaranteed by Theorem 4.1 and Corollary 4.1 of that the Christoffel-Darboux pairing 
between such solutions and any solution of the unbarred DDD system does not depend on x, N ^ or the 
deformation parameters determining the two potentials. It is of interest therefore to compute this pairing 
for the explicit solutions at hand. 

The solutions of the unbarred DDD equations are obtained by taking suitable windows in the wave- 
vector solutions of eq. (2.2.7), consisting of the quasipolynomials vj/'-*'-' (x) and the solutions '^''^\x) 



associated to the f'^^^ (eq. ^.3.16| ). It should be clear that the diagonalization of the pairing depends 



on a careful choice of the contours of integration in eqs. ( 2.3.2 ), ( 2.3.16 ). Here we prove that the choice 



of contours that diagonalizes the pairing ( 2.1.58 ) is linked to the notion of dual steepest descent-ascent 
contours, whose definition is given here and will be needed again in the following section. 
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Definition 2.4 The steepest descent contours (SDC's) and the dual steepest ascent contours (SAC's) for 
integrals of the form 



I{x):= jdye-^^^'^y^-'^y^H^y) , I{x) := Jdy e^'-^^'-y^-'^y^ H{y) 



(2.4.1) 



respectively, passing through any saddle point yk{x), k = . . . ^2 — 1 with H{y) at most of exponential 
type, are the contours 7^ and 7^ , respectively, uniquely defined by 

7fc:=<!yeC; %V2{y)-xy)^^{V2{yk{x))~xyk{x)) ,'R{V2{y)) ^-^^+^ ) (2.4.2) 

7fc:=<(2/eC; ^V^iy) - xy) = ^ (V2{yk{x)) - xyk{x)) ^^{y^iy)) ^-^^-^ ) . (2.4.3) 

Here yk{x) denotes one of the d2 branches of the solution to the algebraic equation 

V^{y) = X , (2.4.4) 

which behaves like 

yk{x) - {vd^+iy^uj''x^ , Lj-.^e^ , (2.4.5) 

X — ^oo 

where x~^ denotes the principal d2-th root of x. (Note that the homology class of the SDC's and SAC's 
becomes a constant for \x\ sufficiently large along a generic ray (see Fig. It will be proved in Section 



3.i that the homology class is also locally constant with respect to the angle of the ray. 
For the following, we also define the sectors 



Sk:^\x^<C, arg(x) G ( ^dy + , , -^y + 7^^^ ]\ , for ^2 odd 

5fc := <j a; e C , arg(a;) G ( --dy + , ~"9y + ''^"j^,^^" ) f- , for 0^2 even 



(2fc- l)7r (2fc + l)7r ' 

2(d2 + l)' ^^^ 2(d2 + l), 
fcTT (fc + l)7r^ 

fc = 0,...,2d2 + l , (2.4.6) 

(where dy was defined in eq. ( ^.1.5| )), and denote by TZk, k — 0,. . .2^2 + 1 the rays separating them, 
counting counterclockwise, starting from Sq. These will be shown in the following section to be inter- 
pretable as Stokes' rays, defining a part of the Stokes' sectors in our Riemann-Hilbert problem. It will 



be shown in Sect. 3.1 that if x approaches 00 along a ray distinct from the 7?.fe's, then the homology class 



of the SDC's and SAC's in Dcf. 2.4 is wcU defined. We then have 



Proposition 2.7 Fix any (non-Stokes) ray {arg(a;) — a — const.} ^ TZu and the corresponding limiting 
(for large \x\) homology class of the dual steepest descent- ascent contours 7^ and ^k, k — l,...d2 of 



Def. 2.4. With this choice of the contours of integration, use the SDC's to redefine the Fourier-Laplace 
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Figure 2: Example of the asymptotic (as |a;| — > oo) SDC's 7 (solid) and SAC's 7 (line-dot-line) in the 
y-plane for a potential of degree d2 + 1 ~ 12 with positive leading coefficient (wda+i > 0) ^^^d for the 
non-Stokes line arg(a;) — 0. The numbers at the vertices of the endecagon label the c?2 = H critical points 
and the pair of dual SDC and SAC passing through each them. 



transform ( 2.3.^ ), and the SAC's to redefine the f^'^^s of eg. ( 2.3.11 ), and the associated wave-vectors 
^C^) modC^^"K Then 



C7"/5 $(")(x)A*^'''(a;) 

N 







Ld2 



d2 + l 1 



N > d2 



up to the addition of a suitable multiple 0/ vp'-^^ (x) to the solutions '^^''\x). 



Proof. 



(2.4.7) 



For P = the statement follows from eq. ( 2.2.36| ) (where the constant c equals \/ho). For a ^ ^ (3 

we use the fact that the solutions ^Ir'-'^^x) associated to /'^'^^ are defined only up to the addition of a 

00 

multiple of v]/^'^)(x). Indeed, we know from the general theory that the pairing of the two solutions is 

00 

constant, and we may use this freedom to normalize the constant to zero by adding a suitable multiple 
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of the homogeneous solution of eq. (2.2.7) \I>^°-'(a;) (i.e., the quasipolynomials), which is "orthogonal" to 



the solutions $^"•'(2;), a 7^ 0, as follows again from eq. (2.2.36). 

00 

Finally fork = a^Q^/3 — j, C''^ equals the bilinear concomitant of the corresponding functions 



(eq. 2.2.34). Let us consider x belonging to a fixed ray, and choose a basis of steepest descent contours 
7fc and steepest ascent contours 7^ (whose homology does not change as x — > cxd on the fixed ray for |a;| 
big enoug h). We use formula ( |2.2.34| ) with 



dye 



-l{V2{y)~-xy) 



'7j V ^0 J'jk 

The respective asymptotic behaviors for large computed by the saddle-point method, are 



/(^")(x) 



2mh 
1 



g- Ji (V-2 (yk (x))-xyk (x)) 



' -27rfi 



(l + 0(A-i)) 



(1 + 0(A-')) , 



(2.4.8) 



(2.4.9) 



(2.4.10) 



where yk{x) are as in Def. 2.4. The bilinear concomitant is a constant, so it must vanish for j ^ k since 



the exponential parts of the asymptotic forms in eqs. ( 2.4.9 ), ( 2.4.10 ) cannot give a nonzero constant 
when multiplied together. 

For j — k, the bilinear concomitant is given by the integral 



1 f - - ,/ ^i(V2(y)-V2(v')-x(y-y')) y-liy) - y2{y') 

y-y' 



2nrh .1^^^^^ 



dy Ay' c'' 



1 



^Ti(V2(yk(x))-xyk{x)) 



2Trh 



2inh \ 



-2T:h 



{l + OiX-'))V^'iy,ix)) 



1 . 



(2.4.11) 



This concludes the proof. Q. E. D. 

— (y) (v) 

Had we chosen the contours as F^ and Fj^ rather than the steepest descent contours, we would 
have had a constant matrix for C"^ . Notice that the pairing of these integrals is also independent of 
the deformation parameters determining Vi, V2, as well as of the choice of the integer N defining the 
window. As shown in this can always be accomplished through a suitable choice of basis. But here 
we have explicitly shown how this occurs for the particular normalizations chosen in the integrals. 



3 Asymptotic behavior at infinity and Riemann— Hilbert prob- 
lem 
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3.1 Stokes sectors and sectorial asymptotics 



Given the duality between the ODE's involving Z?! and Z?^ imphed by eq. (|2.1.60|) , it is clear that the 



Stokes matrices around the irregular singularity at a; = oo for systems of the form ( 2.1.3^ ) and (2.1.39) 
are related. Using the explicit integral representations of the fundamental solutions, we can determine 
the asymptotic behavior from these integral representations by saddle point methods. The solutions for 

N 

the D^i system are simpler to analyze and do not involve the problems of divergent integrals discussed in 



the section p.3| (and the Appendix). But in principle one could consider the asymptotic behavior of the 



solutions to the system (2.1.38) directly, by taking suitable windows in the wave- vector solutions of eq. 



(2.2.7) 



Remark 3.1 We have seen in Sect. 2.4 that one can choose the wave-vector solutions to eqs. (2.2.7), 



(2.2.1 ) so that the Christoffel-Darboux pairing between the dual windows of fundamental solutions implies 
the identity 

N N 

ld,+u N>d2. (3.1.1) 



N N 

*(a;) A*(x) 



N 



N 



Therefore, if we formulate a Riemann-Hilhert problem for $(a;), we can immediately derive the corre- 
sponding one for ^(x). Since the asymptotic forms are mutual inverses, the Stokes (and jump) matrices 



N 



N 



for the one must be the inverses of those for the other. (In our conventions, the Stokes matrices for ^(x) 
act on the left, while for "^{x) they act on the right.) The formulation for the dual wave-functions is 

N 

considerably easier, so this is what we analyze here in detail. 

N 

From the general theory of ODE's, since the matrix I2.iix) is of degree di, one would expect c?i + 1 
Stokes sectors. However, this is true only if the leading term of the matrix has a nondegenerate spectrum. 
In the case at hand, however, we have 



N 

D,{x) 



1 











+ 0{x'^'-^) 



(3.1.2) 



Since the spectrum of the leading term has a d2-fold degeneracy, we have more complicated asymptotic 
behavior and the occurrence of more Stokes sectors. 



A few more preparatory remarks are required. From the discussion in sections 2.2 and 2.3 we obtain 



that the fundamental system for the differential-difference equations specified in Lemma 2.2 is provided 

by 



N 

*(a;) 



$(o)(a;) 



N>di + l, 



(3.1.3) 
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N 



are windows constructed from the wave- vectors defined in Proposition 2.5. Given 



where {^^ 'ix)}k=o.,. 

these integral representations of the solutions, the asymptotic behavior and the Riemann-Hilbert problem 
can be determined by means of the steepest descent method. 

N N 

While {x)}k=i,...d2 are entire functions, the integral defining the window '{x) defines in fact 
a piecewise analytic vector function. Its domains of analyticity are the di + 1 connected components 

(x) 

in which the x-plane is partitioned by the contours Tj , j — 1, . . . rfi. We denote by Vj the connected 

(x) 

domain to the right of the contours for j = 1, . . . , di and by Vq the domain to the left of all contours. 
(See figure for the case di — 7.) The fundamental piecewise-analytic solution defined in eq. 




Figure 3: Example of the two possible choices for the domains of definition of (t)^^"^ (in the x-plane) in the 
case di — 7 



( 3.1.3| ) then satisfies the jump equations 



N 



2i7r>!r'"'i 2i7r>!r'"'2 
1 
1 



2i7r>f^^''2 



1 



N 



(3.1.4) 



where the subscripts +, _ denote the limiting values from the right or the left, respectively, with respect 
to the orientation of the contour. Since we can arbitrarily deform the contours in the finite part of the 
x-plane, we can, by retracting the di contours to the origin arrange that the di + 1 regions all become 

(x) 

wedge-shaped sectors. Using this freedom in the choice of the contours F^- , denote by fi — 0, ... ,di 
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the oriented rays starting from the origin and going to infinity in the sectors defined in (2.1.4). From 
this point on we choose the contours T\ ' as follows 



In this way all the regions, including I?o, become wedge-shaped sectors (see Fig. 
(equivalent) jump discontinuities are 



N 



G^*_(x) 



1 2i7rJp^i 2iTrJfj,2 



,/i = 0, 






■■,di, j 



1, 



22, ^ 



0,j 






1 



(3.1.5) 

The corresponding 



N 

*_(x) 







(3.1.6) 



In order to formulate the complete RH problem, we need to supplement this discontinuity data with 
the sectorial asymptotics around the irregular singularity at x = oo and the Stokes matrices. In doing so 
one should be careful that the lines for which the discontinuities are defined do not coincide with any 



of the Stokes' lines. We can always arrange this by perturbing the rays within the same sector S. 



(x) 
'2fi ■ 



N 



It should be clear that each of the piecewise analytic functions $(a;) can alternatively be analytically 

(x) 

continued to entire functions, since the contours can be deformed arbitrarily in the finite part of the 
x-plane. Therefore the "discontinuities" in the definition of the Hilbert integral are just apparent and 
have an intrinsic meaning only when studying the asymptotic behavior at infinity. Indeed in the final 
formulation of our Riemann-Hilbert problem we prefer to let the lines also play the roles of Stokes' 
lines. 

Proposition 3.1 (Sectorial Asymptotics) In each of the sectors around x = oo with boundaries given 



by the Stokes' lines L^, /i = 0, . . . , di and TZk, k — 0, 



N N N 

fi-r Mx) = Mx) DJx) 
ax 



. 2d2 + 1 defined right after ( 2.4-(\ ) ) the system 

(3.1.7) 



possesses a solution whose leading asymptotic form at x = oo coincides within this sector with the following 
formal asymptotic expansion: 

*'^/orm(a;)~e*^(-)VFa;Gy(a;^) (3.1.8) 

where Y — Yq + 0{x^^^'^^) is a matrix-valued function analytic at infinity, Yq is a diagonal invertible 
matrix (specified in the proof) and 



d2tj 



J + 1 



d2 + l-j , . , 

x^^n'^''+^-^ +Vi{x)E 



(3.1.9) 
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■ 1 







" 











• • • 


W := 









... ^2d. 




. 




^2d, 


• • • w''-" _ 



G := diag ~7V, ■ 



1 _ AT . 

2 2 



-42 iV- 



: 

E 
to 



diag(0,tj,w^ • ■ • 
:diag(l,0,...,0) 



1 



h := - 



1 Vd. 



^2 I'rfa + l 



(3.1.10) 



(3.1.11) 

(3.1.12) 
(3.1.13) 
(3.1.14) 

(3.1.15) 



j — 1 a=oo 

Proof. In any given sector Sk bounded by the lines TZk~i and TZk (eq. 2A.(\ ) we can choose a basis of 
steepest descent contours 7^, k — 1, . . . , 1^2. The reason why the 5fe's are Stokes sectors and the proper 
construction of the steepest descent contours is delayed to the discussion of the Stokes' matrices in Sect. 
3.2. We Fourier-Laplace transform the quasipolynomials 4>n{y) along these contours in order to obtain 
the functions </>^(a;)- Notice that they are not necessarily the same as the previously introduced 0^'^''(x)'s 
since the steepest descent contours do not necessarily coincide with the contours F^^^ defined previously. 
However they are suitable linear combinations with integer coefficients of such (j)^^^ (2;)'s since the choice of 
the steepest descent contours is just a different basis in the homology space of the j/-plane. Now consider 
the asymptotic expansions for 



"P. 



(k) 



(x) := I dye 

Ik 



(3.1.16) 



Here we use the notation ip rather than (p to stress that these are Fourier-Laplace transforms along 
contours of a homology class equivalent to SDC's. The leading asymptotic term in the sector Sk is given 
at the critical point of the exponent V2{y) — xy corresponding to the steepest descent contour ^k- That 
is, we must compute V2{y) ~ xy near a solution to: 



Vi{y) -x = 



(3.1.17) 



asymptotically as a; ^ 00 within the specified sector. Let us solve eq.( |3.1.17 ) in a series expansion in the 
local parameter at 00 given by one determination, A, of the (i2-th root of x: 



Vd.+iy"' + + . . . = Viiy) = x:= X"- 



(3.1.18) 
(3.1.19) 
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We then have the formulae (recalhng A = (V2 (y)) '^^ ) 



.^^dA 1 Vi'(y) ^ 

h = res y(A)— — res y——dy : 

A=oo A d2 y=°o y2(y) 02 ^^^^2+l 



1 



(3.1.20) 



ti = res 



A^' 



J-1 



y'(A)dA = - 



1 



j - 1 y= 



res 



(^2(2/)) dy, J 



As before, denote by {yk{x)}k=o,...d2+i the d2 solutions of the equation V2{y) = x, which are solved 
in( 3.1.20 ) by a series in the d2-th root of x. We then have in a neighborhood oi x = 00 

yk{x')dx' = -d2 / y(A')A"^^-MA' = - ^ ' ^_ Xf^'-' - , (3.1.21) 

J— 

where A^ oj'^A and Cfc is a constant depending only on the coefficients of V2 and the branch of the 
solution yk{x). This formula is proved by taking the derivative (with respect to a;) of both sides and 
using the defining equation for yk{x). Notice that there is no logarithmic contribution since ^^2+1 = 
as follows immediately from the residue formula ( 3.1.20| ). The different saddle points are computed by 
replacing A with Afc := lo'^X. Substituting into the integral representation of the functions ip';^^ (x), we get 

1 



(pW(a;) 

— n ^ ^ 

e ft 
— ' 

Ik. 

e R 



d2/e-*(^=^(^)-="^V„(y) 

Ik 



.*/^!^.(-')d^V„(y(Afe)) / e-iKn'{y(^k))e^^ 



J^2"=(-')d-a„(y(Afe)). 



27rft 



vi'{y{^k)) 



Differentiating V2{y{Xk)) = X^f^ implicitly we obtain the relation 



Vi'{y{Xk)) 



d.X 



2^k 



d2~l 



y'(Xk) 

where y'(A) means differentiation with respect to A. Therefore we obtain 



(3.1.22) 
(3.1.23) 
(3.1.24) 

(3.1.25) 



2TTh , „ _ ''2-1 

-r—Xk ' V, 

hr, 



2T1-1 

2d2 



(oWf^) ~ e-f'-'iV2iVk{x))~xyk(x)) '^njykjx)) / 

^ ^ V ^2" 



2Trh 



iVkix)) 



2-Kh 

d2hn 



Q-h ^iV2{yk{x))-xykix)) ^^n- 



/^e*/^-(^')-a„(y(A.))./M^ 



Ck 



exp 



TAfc 



ft ^ (i2 + 1 - 7 



(3.1.26) 
(3.1.27) 

(l + 0(A-i)) , (3.1.28) 
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where 



Ck ■■= C r^ 



2Trh 



(3.1.29) 



Note that the fuU series (3.1.21) which should appear in the exponent of (3.1.28) has been truncated to 
j < 6,2 because the terms corresponding to j > d2 + 1 contribute to negative powers in the exponential 
and hence give a (l + ©(A^^)) term and, as remarked above, there is no term for j = ^2 + 1 since 
id2+i = 0- 

On the other hand the functions (j)n'^ [x) have the following asymptotic expansion as a: ^ oo within 

0(")(.):=e*-^(^)/d.d,^''""'^"^"^^^^"(^) 



the sectors: 



'=-1 'dsdys'^e-'^^'^i^^V (y) 



fe=0 



(3.1.30) 



JV 



Therefore the matrix of leading terms of #(a:) as defined in eq. (3.1.3) is given by 



diag(l,Ci,---CdJe*^(-) 



7 d.2 —N 

A2^-^-* 



■diag 



_ JV_ 



N + d2 

riN+d2-iVd2+i ^ 
Ai^-^+* 

A2^-^+* 



\d2 



N+d2 
d2 



\/hN- 



Af+d2-l 



(3.1.31) 



The determinant of the Vandermonde-like matrix is very simple: by computing it along the first row 
one realizes that only the first and last minors are not zero. Indeed for all other minors the corresponding 
submatrix has the first and last column proportional. The first minor is a constant in x while the last is 
of order x^'^'^. Therefore we can write 



hN-iVd^i '^^ X' 
Ai 



-N 



N + d2 



N- 



N-i 



1 ^2 



-N~d2 



Ai 



N- 



1 I ti2 



Ad- 



. 1 _L ''S 
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Ai 



Ad, 






1 d2 



Ad; 



1 I •'■2 



Ai 



1 _ d2 



N- 



1 _ ^2 



^^2 



AT- 1 



Ai 



1 I <i2 



TV- 



_ 1 , ti2 



^d2 



1 





[J ^2 2 



2 2 



xdiag (/iAr^ii^d^+i^^""^, A^+5-^, . . . , A^-^+^) (1 + 0{x-^)) 



(3.1.32) 



{l + 0{x-^)) (3.1.33) 



(3.1.34) 



When inserting this into the asymptotic form, we see that, up to factoring the constant invertible (diag- 
onal) matrix on the left 



diag(l, Cic.^+^-*, C2..2(^+i/2-''^/^), • • • , Cd2) 



(3.1.35) 



which is irrelevant for the asymptotics and depends on in a rather trivial manner, we obtain a solution 
with the asymptotic form 



N 1 / 4^ 

^foTm{x) ~ exp - EVi{x) + 2^ 



J=0 



d2 -i + i 



-n 



d2+-l-J 



w 



( N AT M _[_ 1 ^2 

hN-iVd^+i^^x" ,A +2- — 



,A^-^+*') diae 



(l + 0(A-i)) 



diag y/hN-ix 



-AT Wd2+1 A" + 2 2 



exp ( ^WT{x)W-^ ] x-^diag ( ^/h^i, 



IN 

_ N + l 



^^"^ ;; ^ ^ M (i+^(A-)) 

_N+d2 \ 



(1 + C'(A-^)), (3.1.36) 



where W is the matrix defined in eq. (3.1.10). Note that W ^ftW is just the permutation matrix (in 
the subblock). Q. E. D. 
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3.2 Stokes Matrices for the Fourier— Laplace transforms 



The fundamental solution of the system D_i is formed from d2 Fourier-Laplace transforms (a;) , k 



N, 



l,...c?2, and one Hilbert-Fourier-Laplace transform (a 



The asymptotic behavior of the d2 FL 

0,...2d2 + l 



transforms is analyzed by means of the steepest descent method in each of the sectors Sk , k 

N 

( 2.4.6 ), while the behavior of the (piecewise) analytic function '{x) is obtained (in each P^) from eq. 



(y) 



Ay) 



( 3.1.30 ). The computation is achieved by expressing the change of homology basis from the F-j^' 
contours to the steepest descent contours. In order to simplify the analysis of the Stokes matrices we 
point out that there is no essential loss of generality in assuming V2{y) — U'^^'^^ dT+i ' I'^deed, we are 
concerned with just the homology classes of the SDC's, and as a; — > oo the ^2 solutions of the equation 



1^2' (y) = X entering Def. 2.4 are distinct and asymptotic to the d2 roots of a; (up to a nonzero factor). The 
particular choice of the leading coefficient is also essentially irrelevant. If we choose a different coefRcient, 
we must just appropriately rotate hy -d — a.rg{vd2+i) / {d2 + 1) counterclockwise in the pictures to follow, 
but without any essential difference. Therefore, we proceed with Wdg+i set equal to unity. 

With these simplifications, the Stokes phenomenon can be studied directly on the integrals 



dye 



-^{^Tr-^y) 



dz exp 



-A 



d+1 



where 



A:=-N . 



(3.2.1) 



(3.2.2) 



a := arg(a;) , 

and, in order to avoid too many subscripts in the formulae to follow, we have here, and for the remainder 
of this section, have set d = d2- Disregarding the positive factor \x\^ , which is inessential for these 
considerations, the integrals in ( |3.2.2| ) can be written as 



dz exp 



-A 



1 



where z = Z{s) is the D + 1-valued inverse to 



S[z) 



Z{s). 



(3.2.3) 



(3.2.4) 



d+1 

It defines a (d + l)-fold covering of the s-plane branching around the points (z, s) whose projection on 
the s-plane are the d critical values 

d ■ (d + l)a . 



(a) = - 



1 



J = 0, . 



,d-l. 



(3.2.5) 



In realizing this d + 1-fold covering, we take the branch cuts on the s-plane to be the rays 3(s) — 
3(sci'') =const, extending to 5R(s) = -f-oo. As A ^ +00 the integrals ( 3.2.2| ) have leading asymptotic 
behavior that depends only on the critical values of the map s(z) and on the homology class of the 
contour. 
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We now return to the computation of the Stokes' lines. By the definition of the SDC's 7^ (Def. 
2.4, with V2{y) now taken as just y'^^^ /{d + 1)), their image in the s-plane consists of contours which 



come from 5R(s) = +00 on one side of the branch-cut (and on the appropriate sheet) and go back to 
5R(s) = +CX) on the other side of the branch-cut, on the same sheet. (For the SAC's, we choose cuts 
extending to 3fi(s) = —00.) The cuts on the s-plane may overlap only for those values of a = arg(x) for 
which the imaginary parts of two different critical values sir (a) and sii' (a) coincide. A straightforward 
computation, with our simplifying assumptions on V2 (y) , yields the lines separating the sectors Sk to be 
those defined in eqs. ( 2.4.6| ) (with ^2 replaced by d and d — 0). 



We also need the following. 

Definition 3.1 For a given sector S , of width A < tt, centered around a ray arg(x) = ao, the dual 
sector iS^ is the sector centered around the ray arg(a;) = — ao -|- tt, with width n — A. 

(v) 

For a; ^ 00 in each sector the SDC's are constant integral linear combinations of the contours F^ s. When 
X crosses the Stokes line between two adjacent sectors, the homology of the SDC's changes discontinuously. 

(k) 

We denote the SDC's relative to the sector 5^ by 7] , j = . . . d — 1, and denote the column vector with 
these as entries 7^'^-', and similarly, let f denote the column vector with entries {r'f^}j=i...d- Denoting 
the matrix of change of basis by C^, we have 

7W -Cfcf ,Cfc e GL(d,Z) . (3.2.6) 

Our first objective is to compute these matrices Ck- For each fixed generic a (i.e. away from the 
Stokes' lines), we can construct a diagram (essentially a Hurwitz diagram) which describes the sheet 
structure of the inverse map z — Z{s). We draw d + \ identical ordered d-gons each representing a 
copy of the s-plane and whose (labeled) vertices represent the projections of the d critical values s\^) . 
Two vertices with the same label of two different d-gons are joined by a segment if the two sheets are 
glued together along a horizontal branch-cut originating at the corresponding critical value and going to 
3fi(s) = +QO. Since all the branch-points of the inverse map are of order 2, there are at most two sheets 
glued along each cut. Furthermore we give an orientation to the segments (represented by an arrow) with 
the understanding that this gives an orientation to the corresponding SDC. The convention is that an 
arrow going from sheet j to sheet k means that the SDC runs on sheet j coming from 5ft(s) = +00 below 
the cut and goes back above the same cut (or, what is "homologically" the same, the contour runs on 
sheet k coming from -f-oo above the cut and returns to -|-oo below it). 

The diagram can be uniquely associated to a matrix of size d x (d -I- 1), in which each row 
corresponds to a SDC and each column to a sheet. The matrix element {Qk)ij is taken to equal: 

— 1 if the i-th SDC points to the j-th sheet, 
1 if the i-th SDC originates on the j-th sheet. 
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Figure 4: Example of the sheet structure for the case d = 11, a = 0. The contours depicted in the figure 

are the SDC's in the z plane. The second image represents the contours after incrementing a by 7r/6, i.e., 
after crossing two Stokes lines. The labeling of the contours (smaller numbers in the top figures) is by the 
subscripts of the 7fc's and 7fc's passing through the critical points zir^ := e'^o;'^. This also corresponds 
to the numbering of the vertices of the endecagons representing the critical values sif' := —■^^e^^co''. 
The different sheets of the s-plane are mapped onto the connected components in which the 2;-plane is 
cut by the SDC's. These are the images in the 2;-plane of the cuts in the s-plane. The labeling of each 
sheet is given by the boldface numbers and corresponds to the labels inside the endecagons appearing in 
the extended Hurwitz diagrams below the corresponding 2:-planes. The intermediate Hurwitz diagram 
represents the gluing of the sheets after crossing the intermediate Stokes line, i.e., after incrementing a 
by 7r/12. 

otherwise. 

Hence each row of the matrix has exactly one +1 and one —1 entry. 
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In Fig. 1^, we indicate, for the case d = 11, a = 0, the SDS contours and sheet structure represented 
by the extended Hurwitz diagram. The matrix corresponding to this diagram is: 



0123456789 10 11 

1 -1 0000000000' 

1 -1 000000000 

1 -1 00000000 

1-1 

1 -1 

1 0-1 

1 0-1 
0-1 1 
000000 0-1 1 

0000000 0-1 1 

1 00000000 0-1 





1 

2 
3 
4 
5 
6 
7 
8 
9 

10 



(3.2.7) 



where the boldface numbers labehng the columns are the corresponding labels for the sheets, while the 
numbers labeling the rows are the labels of the SDC's. Since a — arg(a;) ranges within a fixed sector 5^, 
the diagram does not change topology and the corresponding matrix Qk remains unchanged. 

We can now describe how a given diagram changes when a ~ arg(a;) crosses the line between two 
adjacent sectors Sk and Sk+i (counterclockwise). These lines correspond precisely to the values of a for 
which two distinct critical values have the same imaginary parts (so that the cuts may overlap if they 
are on the same sheet). We leave to the reader to check that these lines are precisely the boundaries of 



the Stokes sectors Sk defined in eq. (2.4.6). As a increases by 7r/((i + 1) from Sk to iS^+i. the d-gons 
rotate by 7r/d. In this process the connections between the sheets change according to the following rule: 
if the branch-point Pj on sheet r crosses the cut originating from a different branch-point Ph on the 
same sheet (on the left of Pj on sheet r, and hence Pj crosses the cut from below as it moves upwards), 
then the Pj (and its cut) jumps to the sheet s which is glued to sheet r along the cut originating at 
Ph (see Fig. H). Diagrammatically, the tip (or the tail) of the corresponding arrow moves from one 
d-gon to another one connected along the vertex h. In terms of the matrix Qk, the j-th row reflects 
along the hyperplane orthogonal to the /i-th row. The homology class of the SDC 7^'^'' will then change 
because the branch-cut attached to Pj which "emerges" from the branch cut attached to Ph "extracts" a 
contribution proportional to 7^ . The proportionality factor is ±1 depending on the relative orientations. 
The corresponding SDC's 7^'^'' and 7^'^^ are related to the SDC's 7,-'^^"^' and 7^*'''^^^ by the relations 



„(fc+i) 



(fc+i) 



„(fc) 



(k) 



(k) 



(3.2.8) 



ik) (k) 

where the incidence number Chj — ejh is 1 if the SDC's 7^ ,7/j have the opposite orientation and —1 if 
they have the same orientation. Alternatively, the incidence number is just the (standard) inner product 
of the corresponding rows h,j of the matrix Qk- 
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/ Sheet r 


/ 


7 






/ 3 











/ Sheet r 














/ Sheet s 


After crossing 



Before crossing 

Figure 5: The crossing of two critical values and the corresponding cuts. 



We denote by M^. the dx d matrix which expresses the changes in the SDC's described by eq. ( |3.2 
through the relation 

-fC^+i) = Mfc-fW . (3.2.9) 

These will be seen presently to be precisely the Stokes matrices for the passage between these sectors. 
One can check that the matrices Qk and Qk+i are related by 



Qk+i = Mr 



(3.2.10) 



Therefore we can reconstruct the matrices once we have an initial diagram representing the sheet 
structure. 

Before constructing the initial diagram, we note that when a increases by 27r/(d + 1) (i.e. when we 
cross two Stokes lines), the steepest descent (unoriented) contours are like the original ones, but rotated 
by the same amount clockwise. That is, the unoriented diagram (i.e. forgetting the orientation of the 
SDC's) is the same, up to cyclically permuting the labels of the {d + 1) sheets and the d cuts. Note also 
that the critical points rotate by ^^^^^^^ and the critical values by ^ counterclockwise. Indeed we have 



27T 

d+l 



d+1 



e*"e' "i+i z = s(e' ''+^ z; a) 



(3.2.11) 

As for the orientation of the SDC's relative to the new labeling, the d-th SDC passing through zif 
reverses its orientation relative to the (oriented) SDC's obtained by just rotating the initial SDC's (see 
Figure This implies that the matrices Qk representing the diagrams in the various sectors and the 
Stokes' matrices Mk satisfy the recursion relations 



Qk+2=P ^Qk&d+i , Mk+2=p-Mk-p ^ 



(3.2.12) 
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It is therefore only necessary to compute Qq, Qi and Mq, Mi. 

(v) 

Notice also that the wedge contours r!j , j — 0, . . . ,d are in 1 — 1 correspondence with the sheets 
of the map Z{s). Indeed, each connected component to the right of the wedge contours contains only 
one sector 5*2^^1 in the z-plane which corresponds to the sector 5R(s) < in the s-plane. Thus, the same 



argument used to arrive at the recursion ( 3.2.12 ) proves that 



Cfe+2 ^pCkP 



P 



(3.2.13) 
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-1 • 


• • -1 
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e GL(d, Z) , 



where both the matrix &d+i in eq. ( 3.2.12 ) and P in cq. ( 3.2.13 ) generate a representation of the cyclic 
group "Ld+i (although P is of size d x d). The matrix P is just the generator of Z^+i on the "wedge" 
contours; i.e. with the additional constraint X^iLo ^i^^ ~ 

We now describe how to derive the initial diagram - for example, how Figure ^ is obtained. We start 
by noting that the "wedge" contours F^,''^ enclose the sector S2k-i of width 7r/(d + 1). It is to see that 
the corresponding integrals are (more than) exponentially decreasing in the dual sector, since 



.l{y-'+^/id+l)-xy)^y 



< exp(|a;|M) 



dy 



(3.2.14) 



where M is the supremum of -^^di{xy) as y goes along F 



(y) 



Now the contour F^^-* can be deformed so as 



to approach the sector S2^_^ as closely as we wish. Then the constant M is finite and negative if x lies 
within the dual sector S^,^^ . 

We now consider the case of odd d, leaving the easy generalization to even d's to the reader. (The 
only difference is that for odd d, a = arg(a;) = is not a Stokes' line. To study the case even d case, one 
should choose a convenient initial value of a (e.g. a = e^lora = 7r/2(c?+ 1), which corresponds to an 
anti-Stokes line). 

Let us start with the (non-Stokes) value a ~ and focus on any of the SDC's attached to a critical 
value lying in the right s-plane for this value of a (remembering that d is assumed to be odd). Since 
the real part of such a critical value is positive, the corresponding integrals decrease as exp(— ^j^w^ A), 
(A := -^IxIt") i.e., they are (more than) exponentially suppressed on the line x G ]R_|_. As we increase 
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a, the d-gons rotate by ^^ct counterclockwise. It should be clear (from the previous description of the 
change of homology of the SDC's) that the SDC's attached to such critical values do not change homology 
class while they remain in the right half of the s-plane. This is so because there are no branch points to 
their right which can pass through the branch-cut through the given saddle point (which extends to the 
right) as a is increased. Therefore the corresponding integrals are exponentially suppressed as long as a 
ranges in a corresponding sector of width ^y^rj-Ti' — n — j^t^- This is precisely the width of the dual sector 
to a sector of width i.e., the width of any of the -Sj^^j^'s. By careful inspection of the anti-Stokes 
lines for these integrals^, which correspond to the value of a for which the corresponding critical value is 

real and positive, one concludes that they coincide with an appropriate F^,^'' in the left half of the z-plane. 

(v) 

This argument proves that all "wedge" contours Fj^ lying in the left plane are homological to SDC's. 
(These are the SDC's in the left z-plane in Figure ^.) 

We then take the first critical value lying in the left half of the s-plane (in Figure ^, the SDC number 
9). As we increase a so as to move this critical value to the right half of the s-plane, the corresponding 
SDC can acquire a contribution only from the first SDC in the left half-plane (number 8 in our example). 
As a consequence the corresponding integral is exponentially suppressed in a sector of width tt — 2^^^. 
Considering its anti-Stokes line and its linear independence from the previously identified SDC's, we 
conclude that it must enclose two odd-numbered sectors >S'2fe+i and S2k+3 (in our example this is contour 
9). Proceeding this way we can easily identify the homology classes of all SDC's for a = 0. 

The labeling of the sheets is largely arbitrary. (The choice we have made in the fig. 4 is just for 

(v) (k) 

"aesthetic" reasons). The fixed basis of contours FJ , i — 1, . . .d and the bases 7j , j = 0, . . . d — 1 fc = 
0, ... 2d -I- 1 are related by 

f^"-) = Ckf = A4_i • • • MoCof (3.2.15) 

7(0) = Cof . (3.2.16) 

Therefore only the matrices Mk and the first change of basis matrix Cq are needed. The matrix Co 
can easily be constructed from the initial sheet structure. For our present purposes it is not actually 
necessary to have its general form, since the objects of primary interest are the Stokes matrices, which 
will be seen to be just the M^'s. But to illustrate by example the form that Cq takes, within the basis 



®The anti-Stokes' line for exponential integrals of the type of Def. 2.4 are the lines along which the integrals are most 
rapidly decreasing as |x| ^ 00. 
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we have chosen, for the case d = 11, it is: 



11111100000 
1 1 1 1 
00110000000 
00010000000 
00001000000 
Cn:= 00000100000 (3.2.17) 
00000010000 
00000001000 
00000000100 
00000000110 
1 1 1 1 

It is not difficult to give an exphcit description of the matrices Qq, Qi and Mo, Mi but, for the sake of 
brevity, we will not give it here. It consists of a lengthy but straightforward calculation, which leads 
to the matrices in Table |l| (listed for the cases cZ = 2, . . . , 11), where one can clearly extrapolate to the 
correct rule. 

We now turn our attention to the Stokes matrices. Consider the subblock of the fundamental system 
corresponding to solutions of the D.i ODE given by the Fourier-Laplace transforms. If we denote by 
Y{x) the d2 X ^2 such block whose rows are the integrals of fiy) on the contours and by Yk{x) the 
analogous matrix obtained by integrating over the SDC's, we have 



Yk{x) = CkY{x). 



(3.2.18) 



The Stokes matrices are then given by 



Sk ■= Yk+iYk 



Ck+iCk 



Mk 



(3.2.19) 



These are just the matrices expressing the relative change of homology basis of the SDC's corresponding 
to two consecutive Stokes' sectors. 

In order to complete the description of the Riemann-Hilbert problem, we need to also consider the 
extra solution given by a Hilbert-Fourier-Laplace transform. In doing so, we extend the previously 
computed Stokes matrices Mk to the full fundamental system of d2 + 1 solutions by means of 



1 








Mk 



Mk := 



Summarizing the whole discussion, we have proved the following theorem. 



(3.2.20) 



N 



Theorem 3.1 (Riemann— Hilbert Problem) There exists a fundamental system of solutions, J&(x) 
which is analytic where defined, with the properties: 
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Mo = 
Mq = 

Mo = 
Mq = 

Afo = 
Mo = 



Mo = 



Mo = 



Mo = 



Ml = 

Ml = 

Ml = 

Ml = 

Ml = 

Ml = 



Ml = 



Ml = 



Ml = 



Table 1: The first two Stokes matrices for degrees d2 = 2 ... 11. The remaining ones are obtained using 
eq. ( pl^ ). 
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1. In each of the sectors around x = oo defined by removing the lines L^, fi — 0, . . . di and TZk, 

N 

fc = 0, . . . , 2d2 + 1, ^{x) exists, is analytic and invertible, and it can be normalized on the left by a 

N 



constant matrix to have the same asymptotic behavior as ^formix) defined in Proposition 3.1 



2. Crossing the lines L^, the corresponding Stokes matrix is given by as defined in eq. (S.l.t) 



3. Crossing the lines TZk the corresponding Stokes matrix is given by in eq. (3.2.2C), constructed 



according to the algorithm described in section 3. 




Figure 6: Example of the structure of Stokes sectors and discontinuity lines (which can equivalently 
be viewed as Stokes lines) in the x-plane for the case di — 7 and d2 — 11 and both leading coefficients of 
the potentials real and positive. 

For an example of all Stokes' sectors and lines see Fig. ^ and for examples of explicit Stokes' matrices 
Mq , Ml from which all the others are computed easily, see Table 0. 



We conclude this section with the remark that in the formulation of Theorem 3.1 the lines were 



viewed as Stokes' lines. However, we could alternatively have formulated an equivalent Riemann-Hilbert 



problem, in which the lines define discontinuities, according to eq. (3.1.q), by retaining a fixed choice 



of the contours of integration along the boundaries of the wedge sectors, thereby giving rise to genuine 
jump discontinuities across the boundaries. 
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4 Summary and comments on large A'^ asymptotics in multi- 
matrix models 



A complete formulation of the Riemann-Hilbert problem characterizing fundamental systems of solutions 
to the differential - recursion equations satisfied by biorthogonal polynomials associated to 2-matrix 



models with polynomial potentials is provided in Theorem 3.1. The approach derived here can also be 
extended quite straightforwardly to the case of biorthogonal polynomials associated to a finite chain of 
coupled matrices with polynomial potentials, following the lines indicated in the appendix of [Q. 

The Riemann-Hilbert data, consisting essentially of the Stokes matrices at oo, are independent of 
both the integer parameter TV corresponding to the matrix size and the deformation parameters deter- 
mining the potentials; that is, the fundamental solutions constructed are solutions simultaneously of a 
differential-difference generalized isomonodromic deformation problem. This provides the first main step 
towards a rigorous analysis of the N ^ oo, hN = 0{1) limit of the partition function and the Fredholm 
kernels determining the spectral statistics of coupled random matrices (essential, e.g., to the question of 
universality in 2-matrix models) . Such an analysis should follow similar lines to those previously success- 
fully applied to ordinary orthogonal polynomials in the 1-matrix case [|[ |9[ 10, The main 



difference in the 2- (or more) matrix case is that in the double-scaling limit the functional dependence 
of the free energy on the eigenvalue distributions is not as explicit as in the 1-matrix models j2^, p^ . 
It is also clear that the hyperelliptic spectral curve that arises in the solution of the one-matrix model 
must be replaced by a more general algebraic curve, which arises naturally in the spectral duality of the 



spectral curves of |g] (cf. |14| ) 



In order determine the large N asymptotics with the help of the data defining the Riemann-Hilbert 
problem, one should begin with an ansatz that can be checked a posteriori against the given case. In the 
1-matrix case ||l^, |l^ , this was provided by means of hyperelliptic 9-functions. The physical heuristics 
and the basic tools for generating such ansatz were also given in |l^, |2^, ^ . Much of of these heuristics 
can be extended to the 2-matrix case Q, and this will be the subject of a subsequent work [Q. 

A Appendix: Convergent integral representations of 

oo 

In this appendix, we indicate how to overcome the problem in the formal definition of eq. ( |2.3.15D , 
illustrated through an example. The idea follows that used in |^ for the case of cubic potentials^. 

As discussed in section v|/'''^(a;), a natural approach to finding solutions of eq. ( |2.2.7D would be to 

oo 

take the inverse Fourier-Laplace transform of the function 

r — — {Vi{x)—xt) 
Jf(")(y) ■-ei^^'^y^ dxdt— **(a;) , (A-1) 

oo J^r y 



^We wish to thank A. Kapaev for pointing out the problem of divergent integrals appearing in an earlier version of this 
paper, and for helpful discussions. The approach outlined here, for potentials of arbitrary degrees, will be the subject of 
further work, and the results detailed in a subsequent, joint publication 
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along the anti- wedge contours T^^^ ■ Eq. (A-1) defines a piecewise analytic function in the domains 

d2 



k 

k=l y=0 



c.\Urr' = U^- (A-2) 



with jumps across the F^^'' given by 

di 

4-(0)(y+)-^.(0)(y_) = 2z^^>^^''=*^-'")(y) , v&T^^\ (A-3) 

oo oo . ^ oc 

However, either retaining this definition or modifying it through analytic continuation leads to a divergent 
integral when the Fourier-Laplace transform is applied. In order to resolve this difficulty, the idea is to 
consider a suitable contour approaching y = oo within the odd-numbered sectors 5^^*^^^^ and define a 
locally analytic function in a neighborhood of this contour so as to render the Fourier-Laplace transform 
convergent. 

A general approach resolving this problem is not needed for the results presented in this paper since, 
as pointed out, the Riemann-Hilbert problem for '^{x) is more simply obtained by means of duality 

N 



(Remark 3.1). Therefore, we just give an indication of how to proceed by means of an example. 



A.l An Example: di = 7, ^2 = 4 

The example we consider consists of two potentials Vi{x) and V2{y) of degrees 8 and 5 respectively, both 
with positive real leading coefficients. We refer to Fig. |^ as the "the figure" throughout this section. The 

(k) (j) 

Sy sectors are the light-shaded ones while the Sx are the darker-shaded ones. Fix an asymptotic ray 
£c with arg(y) = Q!o(= ^) within an odd-numbered sector, in our example Sy^\ and an e sector around 
it. Let C be the dual sector, i.e., the sector centered around tt — ao(= 0) with width tt — e. Any of the 
functions ^''•'\y) whose wedge contours rj^"* lie in the sector C decay faster than any exponential e^^^^ 

oo 

on the given ray £c- These are the sectors ^5^^'"^'^^'^'^) jj^ Fig. |^. Fix another ray ifi within a different 
odd-numbered sector Sy'^'^'^^'^ in such a way that the dual sector to the e sector around £7^ (denoted by 
TZ) contains all but one of the remaining sectors Sx^'^^'^K In the choice shown in the figure i-ji lies in 
iSg^'' and the sector TZ contains the sectors si'^'^'^^^ Finally, fix a contour F which goes off to infinity 
asymptotic to the two rays £c, (as in figure). 

With these choices, redefine the duality pairing according to the following rule: 

(i) In the x-plane use as basis only the contours F^^'' which fall within one or the other of the sectors 

fa;) 

£, TZ. In our example this leaves out contour F3 , which has not been drawn in the figure. 

(ii) In the y-plane use a basis in which the contours F intersect only one of the F^^''s. In our example this 

(y\ (y\ 

would leave out either Fg or Fq and we have chosen the latter. 
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Figure 7: The figure illustrating the definition of 0(-^^/*'^*ff''*)(y). The contour F must approach oo 

oo 

within the sectors ^j^^j^'s (lighter shaded) and asymptotic to the two rays H-cAn- The terms added to 
in the definition of Q'^^'^f*'' (e(-R*9''*)) must correspond to sectors S^f,^^ which faU in the dual sector 

oo oo oo 

C [TV) indicated by the dashed (dash-dotted) arc. 



Note that we can accomplish an equivalent duality pairing using such bases by exploiting the homo- 
logical equations 

Erf = Eri^^ = o. (A-i-i) 

j=0 fc=0 

Indeed, using the homological equation ( |A-l-l| ), one can avoid one contour and redefine the x's: 

>iT= ><'''^r|") X F*^' = ^''"^f-^^f, (A-1-2) 
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where x^'^ in this example would be 

^i,3 .= _ ^ i ^ 3,0; 3^°'^' := -x^^' . (A-1-3) 

Now define the wave- vector 

0i<=/*(y) :=*(")-*(j/)-2i7r ^ 5^*'3^W,t(y) =; _ 2i7r ^ 3^''3 , (A-1-4) 



1=0,1,2,7 



for y belonging to the part of T inside contour T^^\ (Note the transposition of row vectors to column 
vectors.) In eq. ( A-1-4 ) in general Ic would be the set of indices j corresponding to the contours F^^-* 
contained in sector C. Note that this function decays on the ray ^c because is in the dual sector of 
all the \I> entering eq. (A-1-4). Moreover the piecewise analytic function vj/^"'' decays exponentially in 



each sector S, 



(v) 



As we cross the contour Fg along F, by virtue of the jump condition ( |A-3D , the function is 

oo 

analytically continued to the function 

0«'ff''*(y) := ¥°^'\y) + 2i7r^S^'^3^(^),t(y) _ *(°)'*(2/) + 2i7r^ :^<^i^''{y) , (A-1-5) 

°° j^L °° °° j=4,5,6 °° 

which, due to the choices made above, also decays along the asymptotic direction i-ji. The two functions 
QLef t, Right ^y-^ define together a function @^{y) which is (locally) analytic in a neighborhood of the 

oc oc 

contour F and decays sufficiently fast in either direction to allow us to take an inverse Fourier-Laplace 
transform using 



¥^H^) /dye-^0(r)(2/) 



(A-1-6) 



This wave-vector is a solution of cqs. (2.2.8), (2.2.40) and (2.2.41) as we now show. The equations 

(A-1-7) 
(A-1-8) 
(A-1-9) 

(A-1-10) 



y0(f)(y) ^P0(f)(y)+ y^e*^=(^),0,. 

OO oo L 

hdy &^^\y) = Q&^^\y) + W2{y) [v^e*^^(«), 0, , 

OO oo 

hd^^@(^\y)^V'<@^^\y) 



hd,,@^^\y) = -v'e'^^\y)- 



1 P ' - y' 



J P-y 

are satisfied as a consequence of the similar equations for ^'^\y) and ^^\y), fc = 1, . . . di. From these, 

oo oo 

eqs. ( ^.2.8| ), ( ^.2.4C ) and (2.2.41) for vj/'-'"'' (a;) follow from standard manipulations of the inverse Fourier- 

oo 

Laplace transform which are now made rigorous by the convergence of the integral. In particular this 

dye* ^^^(2')-^^) . (A-l-U) 



solution of (2.2.8) is associated with the function 
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One should repeat the scheme outhned here for other contours F as weh, until they span the same ho- 
mology space spanned by the antiwedge contours {f'-^''}, so as to obtain a basis of solutions to eq. (2.2.7) 



Just to complete the given example we indicate how to choose the other contours and corresponding bases: 

1. The contour asymptotic to the rays arg(?/) = — tt and arg(y) = tt/IO + e. The basis on which to 

j 



re-define the coefficients x (which we now denote by >f) is given by T^^ X F^''^ i ^ 6, j 7^ f and 



the wave vectors are given by 

0i<=/*(y) :=^(0)^*(y)_2i7r ^ J^^3^W^*(y) (A-1-12) 
@ffis/«t(y) := 4^(o)'*(y) -f 2^7^^ k^-'^ ^^^'^'\y) ■ (A-1-13) 

j=3,4,5 

2. The contour asymptotic to the rays arg(y) = 7r/2 + e and arg(j/) = 37r/2 — e. The basis on which 
to re-define the coefficients x is given by v'f'^ x F^-^'', i ^ i, j ^ A. 

0^«/*(y) ¥°^^\y) - 2i7r^ (A-1-14) 
0ffisftt(j^) := 5^(0)'* (y) 2i7r ^ S^'^^ ^^^^t^^) _ (A-1-15) 

j=5,6,7,0 

3. The contour asymptotic to the rays arg(y) = tt/2 + e and arg(j/) = ITtt/IO + e. The basis on which 
to re-define the coefficients x is given by v'f^ x F^-^'', i ^ Q, j ^ Q. 

0^«/*(y) «^W'*(y) - 2i7r ^ S^^'^ ^^W>*(2/) (A-1-16) 

^=1,2,3,4 °° 

QR-ght^y^ :=_5-(o)^*(y) + 2i7r^ JJ^-^ ^^(^■)^*(y) . (A-1-17) 
We leave the details of the general case to subsequent work. 
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